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We consider the interplay of thermal fluctuations and shear on the surface of the domains in 
various systems coarsening under an imposed shear flow. These include systems with nonconserved 
and conserved dynamics, and a conserved order parameter advected by a fluid whose velocity field 
satisfies the Navier-Stokes equation. In each case the equation of motion for the interface height 
reduces to an anisotropic Burgers equation. The scaling exponents that describe the growth and 
coarsening of the interface are calculated exactly in any dimension in the case of conserved and non- 
conserved dynamics. For a fluid-advected conserved order parameter we determine the exponents, 
but we are unable to build a consistent perturbative expansion to support their validity. 



PACS numbers: 05.70.Ln, 05.40. +j, 02.50.-r, 81.10.Aj 



I. INTRODUCTION 

This paper deals with the influence of shear on inter- 
facial fluctuations in phase-ordering or phase-separating 
systems. The primary motivation is the need to under- 
stand the influence of thermal fluctuations on coarsen- 
ing under shear. Thermal fluctuations are not normally 
thought to be important for coarsening systems, as the 
dynamics is controlled by a "strong coupling" , i.e. zero- 
temperature, fixed point and temperature is formally 
an irrelevant perturbation Q|. Under an externally im- 
posed shear flow, however, the growing domains become 
stretched in the flow direction and there is evidence, 
especially in two spatial dimensions, that growth in the 
transverse direction is strongly suppressed This 
raises the possibility that thermal roughening of the in- 
terface might destroy the coarsening state. On the other 
hand, the thermal roughening is itself suppressed by the 
shear flow, so the question of the survival of the coars- 
ening regime to late times rests on a delicate balance 
between these two effects. 

A second motivation for this study emerges from the 
mathematical description of the interfacial fluctuations, 
which takes the form an anisotropic Burgers equation 
P,fio[ . The structure of the equation, and the form of the 
noise correlator, are such that, in a renormalization group 
(RG) analysis, some parameters of the theory are not 
perturbatively renormalized. As a result, certain combi- 
nations of scaling exponents can be determined exactly. 
Remarkably, the number of such combinations is in every 
case equal to the number of unknown exponents, so that 
all scaling exponents can be determined exactly for any 
spatial dimensionality d. 

The structure of the interface equation is very simple. 
If /i(x, t) is the interfacial height relative to the mean 
height, where x is a (d— l)-dimensional vector specifying 
position in the plane parallel to the (mean) interface, and 
t is the time, the equation takes the simple form 



where 7 is the shear rate, and we have taken the shear 
flow to be in the x direction. The linear operator C is 
diagonal in Fourier space, and its eigenvalues A(k) have 
the limiting small-fc form 



A(k) ~ |k| 1+ ^, 







(2) 



d t h + ^hd x h — Ch + r/(x., t) , 



(1) 



In equation ([!]) we have retained only the leading- 
order nonlincarity, which is associated with the shear. 
In this limit, the noise correlator has the same form 
as in the zero-shear case, namely (in Fourier space) 
(?7(k,f)?/(k,i')} - Ik^-^k+k')^-*'), where this par- 
ticular form follows, via the fluctuation-dissipation the- 
orem, from the zero-shear stationary state, P[h(x)] oc 
exp[— const / d d x(\7h) 2 ]. The parameter fx specifies the 
particular dynamical model under consideration. Partic- 
ular cases of physical relevance are fi = 1 (a nonconserved 
order parameter, or 'model A' in the classification of 
Halperin and Hohenberg [0), (1 = 2 (a conserved order 
parameter obeying the Cahn-Hilliard equation, or 'model 
B'), and /1 = (a conserved order parameter coupled to 
hydrodynamic flow in the viscous regime, or 'model H'). 

The derivation and RG analysis of equation (|l|) will 
form the main part of this work. Since the system is 
anisotropic due to the shear, we write x = (x, xj_), where 
x is the coordinate along the flow direction, and xj_ is 
a (d — 2)-dimensional vector perpendicular to the flow. 
There are, in general, three scaling exponents, \i C an< i 
z, defined by the condition that the simultaneous scale 
transformations x — > bx, x.± — > 6^xj_, h — > b x h, and 
t — > b z t leave the interfacial dynamics scale invariant. 
All three will be determined exactly for all physical val- 
ues of \i and for all d. 

The remainder of the paper will consist of a more de- 
tailed discussion of the physical motivation for these cal- 
culations, and the analysis and interpretation of the re- 
sults. The interface equations are derived in section II for 
models A, B, and H. Section III contains the RG anal- 
ysis, while in section IV we discuss the implications of 
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our results for coarsening systems under shear. Section 
V concludes with a summary of our results. 



II. THE INTERFACE EQUATION 

In each case we will start from the relevant Ginzburg- 
Landau equation for the order parameter <p(r, t), and de- 
rive the interface equation by projecting the full equation 
of motion onto the interface. We assume a coarse-grained 
free-energy functional of the Ginzburg-Landau form, 



F[4>] = J d d 7 



(v^) 2 + v(4>) 



(3) 



where V(cj>) is a symmetric double-well potential with 
minima at <p — ilj representing the two equilibrium 
phases. 

For pedagogical purposes we begin with the simplest 
case of the time-dependent Ginzburg-Landau equation 
(or 'model A') which describes phase-ordering in a system 
with a nonconserved scalar order parameter, i.e. Ising-like 
systems such as a twisted-nematic liquid crystal. 



A. Model A 

We will consider a uniform shear flow in the re- 
direction, with the velocity gradient in the y-direction, 
v = 7J/e x , where 7 is the shear strength. The dynamics 
of the the system are governed by the Langevin equation 



d<j) d<j> SF „ . 



V 2 0-y'(0)+£(r,i) , 



(4) 



where the second term on the left-hand side is just v ■ V</>, 
and represents the advection of the order parameter by 
the shear flow. In equation (Q), a kinetic coefficient has 
been absorbed into the timescale, V'{4>) = dV/dcj), and 
£(r, t) is Gaussian white noise with mean zero and corre- 
lator 



{Z(r,t)ar',t'))=2Dd(r-r>)5(t-t'), 



(5) 



where the noise strength D is proportional to the tem- 
perature. 

We now construct an equation for an interface, parallel 
to the flow direction and normal to the velocity gradient, 
separating the equilibrium phases. We are interested in 
the limit where the interface is almost planar, such that 
(Vft.) 2 is typically small, i.e. we are going to systemati- 
cally neglect terms which are smaller by powers of (Vft.) 2 
than the terms we retain. In this limit, the order param- 
eter profile is well represented by the simple form 



where we have written r = (x, y). The function f(u) 
is essentially a step function, with a width given by the 
interfacial width, £o- Its derivative, f'{u), is therefore a 
smeared delta function, which peaks on the interface and 
has width £q- It will be used below as a projector onto 
the interface. 

Substituting equation (^) into equation (Q) gives, with 
u = y-h, 



(dth)f(u) 



-l(u + h)(d x h)f'(u)- 
-(V 2 h)f'(u) - V{f) 
-£[x, u + h(x, t),t] . 



[1 + (vm 2 ]/» 



(7) 



Finally we multiply through by f'(u) and integrate over 
u. Formally we take the integral from —00 to 00, but in 
practice the integral is concentrated in the neighborhood 
of u = 0. Since f"{u)f{u) and V'(f)f'(u) are perfect 
derivatives, these terms drop out. Also the term involv- 
ing u[f'(u)] 2 vanishes by symmetry under the integral. 
The final result, therefore, is 



d t h + jhd x h = V 2 h + T]{x, t). 
The noise term is given by 

T]{x,t) = -(1/(7) / duf'(u)S[x,u + h(x,t),t] 



{r,t)=f[y-h{x,t)] 



(6) 



(8) 



(9) 



where a = j du[f'{u)] 2 is the surface tension. Clearly 
the mean of rj is zero, while use of equation (||) gives its 
correlator as 

(tj(x, t)r,{x', f )> = (2D/a)5(x - *)8(t - t') . (10) 

In the zero-shear limit, 7 = 0, equation (|^) reduces to 
the Edwards- Wilkinson model and has a simple in- 
terpretation. The interfacial free energy functional, to 
lowest order in (V/i) 2 , is F mt = {1/2] )J d d ~ 1 x (Vh) 2 . 
The dynamics (|S|) corresponds to the Langevin equa- 
tion dth = —8Fi n t/8h + 7]. The noise strength 2D /a 
in (^) guarantees the correct stationary distribution, 
P[h] <xexp{-aF[h]/D). 

Before moving on to model B, it is worth noting that 
for the case of zero shear and zero noise the equation 
reduces to simple relaxation. In Fourier space, one has 
dh(k,t) = —k 2 h(k,t), i.e. fluctuations on a length scale 
L ~ 1/k relax on a timescale r{L) ~ L 2 . For a coarsen- 
ing system containing many interfaces, this relation gives 
the timescale, L 2 , for a feature at scale L to relax away, 
and suggests the relation L{t) ~ t 1 / 2 for the coarsening 
length scale, or 'domain scale' in a phase-ordering system. 
This approach to determining coarsening exponents from 
interfacial relaxation rates has been used before 
and the predictions agree with the results obtained from 
other methods M. Indeed, the result is more general 
p"5| . In any system where coarsening proceeds by relax- 
ation of extended defect structures (domain walls, vor- 
tex lines, etc.) the dynamical exponent 2, in the relation 
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Lit) ~ t 1 /* for the coarsening dynamics, is the same as 
that obtained from the relaxation rate, A(k) ~ |k| z , of a 
single defect with a sinusoidal modulation at wavevector 
k. The same general structure will be apparent in the 
study of models B and H. 

B. Model B 

For conserved dynamics, the time-dependent 
Ginzburg-Landau equation is replaced by the Cahn- 
Hilliard-Cook equation (i.e. the noisy Cahn-Hilliard 
equation) which, in the presence of a uniform shear 
flow, reads 

d<j> d<t> -6F . 

= -V 2 [V 2 0-U'(0)]+£(r,i) , (11) 

where a transport coefficient has been absorbed into the 
timescale. The noise correlator is 

(e(r,t)e(r',t')) = -2 J DV 2 5(r-r')5(t-0 . (12) 

As a prelude to further analysis it is convenient to first 
operate on both sides of the equation with the inverse 
of the Laplacian operator (whose meaning will become 
clear below). Making the same long- wavelength approx- 
imation @ as in the treatment of model A gives 

(~W 2 )-\d t h + 1 (u + h)d x h)f'(u) 

= -[1 + (Vh) 2 ]f"(u) + (V 2 h)f(u) - V'{f) 

+ (-V 2 )- 1 e[x, M + / l (x,i),i] . (13) 

Multiplying through by /'(w), and integrating over u as 
before, gives 



duf [u^-V 2 )- 1 f (u)[d t h + 7 (« + h)d x h] 

= aV 2 h + fj(x, t) , (14) 

where the noise is given by 
fj(-x,t) = -[ d M /'H(-V 2 )- 1 ^[x, U + / l (x,t),i] . (15) 

The meaning of the operator (—V 2 ) -1 is as follows. In 
Fourier space one has (— V 2 ) -1 — » (k 2 + q 2 )^ 1 , where 
(k, q) is the vector conjugate to (x, y). Defining, for 
a general function F, G(x, y) = (— V 2 ) _1 F(x, y), its 
Fourier transform, in the (d — l)-dimensional subspace 
spanned by x, is given by 



G{Ky) = 2\k\J d v' e M-m\y~y'\)F(Ky') ■ (16) 

We now use this result to evaluate the left side of equa- 
tion (|T4|). The leading order non-linearity (in h) is given 



by the shear term, so elsewhere in equation (|IJ) we ne- 
glect the distinction between u and y = u + h. It can be 
shown that the leading-order terms omitted in this ap- 
proach are of order h(d x h) 2 . Denoting, for brevity, the 
Fourier transform with respect to x by a subscript k, the 



Fourier transform of the left-side of (14) becomes 



J du J dv exp(-|k| \u - v\)f'(u)f'(v) 



d t hk + ijk x ( v + -^[h 2 ]k 



(17) 



Recalling that f'(u) acts like a delta function at u = 
(of strength 2, which is the discontinuity of the order 
parameter across the interface) equation (|l4| ) simplifies 
to 



d t h k + -jk x [h\ 



-||k| 3 / lk + l|k|ry k (t). (18) 



Consider once more the case of zero shear and zero 
noise. Then equation ( |l8| ) represents simple relaxation, 
with fluctuations on length scale L ~ 1/k relaxing at a 
rate fc 3 , i.e. as k z with z — 3. This is again consistent 
with the known coarsening growth law, L(t) ~ t 1 / 3 , for 
model B §. 

The form of the noise correlator can be extracted from 
equation <^). Using the same simplifications as before 
yields, in Fourier space, 



AD 

(fj k (t)n- k >(t')) = —6 K v6(t-t') 



(19) 



Equation (jig ) has, in real space, precisely the form 
of equation (|l|)^ where the operator C has the small- |k| 
spectrum A(k) ~ |k| 3 , i.e. it has the form (g) with fx = 2. 
Defining r)u.(t) — ^\k.\fjk(t), one recovers equation (|l]) ex- 
actly, with noise correlator 



(rik{t)v-h>V)) = D\k\ <J k , k - S(t - t') 



(20) 



For model A, equation (||) also has the form (|J), but 
with n = 1 in (Q). This suggests that both models be 
viewed as members of a more general class, defined by 
equations ([!]) and (0) with /i general. As discussed in 
the Introduction, the requirement that the equilibrium 
distribution P[h] oc exp[— ^ ^ k k 2 h\Ji_-k\ be recovered 
for 7 = forces the noise correlator to have the form 
fak(t)»7-k'(t')) ~ Ikr" 1 <y k)k / 5{t - t 1 ). Our results © 
and (pC|), for models A and B respectively, satisfy this 
requirement. 
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C. Model H 

The general results relating the form of the spectrum 
(||) of the operator C in (Q) to the exponent z for coarsen- 
ing {Lit) ~ t x l z ), and the form of the noise to the require- 
ment of recovering the correct equilibrium state in zero 
shear, suggests a simple form for the equation of motion 
for an interface in a phase-separating binary fluid in the 
'viscous hydrodynamic' regime. This is the regime de- 
scribed by 'model H' of the Hohenberg-Halperin scheme 
jnj. In this regime, it is known that coarsening proceeds 
linearly in time, Lit) ~ t, corresponding to z = 1 |jl6| . 
This suggests that the interfacial relaxation spectrum is 
given by A(k) ~ |k| for k — > 0, i.e. p = in (||), a re- 
sult which has been confirmed by Shinozaki |L4j . This in 
turn suggests that the interfacial noise correlator should 
have the small k form corresponding to /i = 0, namely 
(77k(t)r?- k '(*')) =D\k\- 1 6 icM ,6(t-t'). 

We now show that these expectations, based on gen- 
eral considerations, are indeed borne out in practice. In 
the absence of thermal noise, the equation of motion for 
the order parameter field takes the form 



at 



(21) 



where fi — 8F/5<f> is the chemical potential and T is a 
transport coefficient. The velocity, v, of the fluid, as- 
sumed incompressible, satisfies the Navier- Stokes equa- 
tion 



M^ + (W)v 



7iV 2 v - Vp - </>V^ 



(22) 



where p and r\ are the density and viscosity of the fluid 
respectively, and p is the pressure. The final term in ( p^ ) 
contains the feedback between the order parameter and 
the fluid velocity. 

The coarsening dynamics of this system is known to 
exhibit three regimes |l6| , |T7| ]: (i) an early time 'diffusive' 
regime, where the hydrodynamics is irrelevant (the fluid 
velocity is much smaller than the typical interface veloc- 
ity) and the model reverts to model B, with coarsening 
scale L(t) ~ t 1 / 3 ; (ii) an intermediate time 'viscous hy- 
drodynamic' regime, where the 'inertial terms' on the left 
side of equation (|2^ ) can be neglected, with Lit) ~ t; 
(iii) a late time 'inertial hydrodynamic' regime where 
the inertial terms dominate the viscous term, 77V 2 v, and 
L(t) ~ t 2 ' 3 . 

Here we focus on the viscous hydrodynamic regime, 
where we can set the left side of (22) to zero. This de- 
fines model H Q|l]]. For simplicity, we will ignore the 
imposed shear flow in the first instance. The pressure 
can be eliminated by using the incompressibility condi- 
tion, V • v = 0, to express the velocity in terms of (j>Vfi. 
Putting the result into (^T|), and adding a noise term 



gives the final equation for model H. Since we are inter- 
ested in the regime where diffusion is negligible, we drop 
the term TV 2 /i to obtain 



dt 



dr' cW(r) T ab (r - r') cW(r') /i(r') + £(r, t) , 

(23) 



where \i — 8F/54> — V'(<p) — V 2 </>, and T ab is the Oseen 
tensor, with Fourier transform 



T ab (k) 



1 



Sab 



rjk 2 V k 2 



(24) 



In equation (P3f) , repeated indices are summed over. The 
form of the noise correlator is dictated by the fluctuation- 
dissipation theorem: 

(f (r, m{r',t')) = 2Dd a <j ) {v) T ab (r - r') 3 b 0(r') 

Xfji{r')S{t-t') , (25) 

where D is the temperature. 

To determine the interface equation we insert the form 
(||) into ( p3] ) to obtain, analogous to (0) 



(d t h)f(u) = J dr' cW(r) T ab {r - r') ^(r 7 ) 

x{(V 2 h)f'{v) + V'[f(v)} 

-[l + {Vh) 2 ]f"(v)} 

-£(x,u + h(x,t),t) , (26) 

where u = y — h(x,t) and v = y' — /i(x', t). It is impor- 
tant to note that the Oseen tensor in real space is only 
defined for d > 2. Therefore, all the following equations 
for model H are only valid for d > 2. 

As in models A and B, the leading term for small h 
comes from the V 2 /i term in the braces. To linear order, 
therefore, we can use a 'flat interface approximation' in 
the terms outside the braces. This means we can write 
V</>(r) = f'{u)e y , V</>(r') = f'(v)e y , where e y is a unit 
vector in the y direction, and T yy becomes the only rele- 
vant element of the Oseen tensor. Multiplying both sides 
of ( p6| ) by f'{u), and integrating over u, yields, to leading 
order in h, 

d t h{x)=(7 J dx'T TO (x-x',0) V 2 h{x) + noise , (27) 

where the integral is over the {d — l)-dimcnsional plane 
of the mean interface. Fourier transforming this result, 
using (p||), gives 



dhk 
dt 



At) 



%(*) > 



(28) 



where k is now a (d— l)-dimensional vector, and we recall 
that d > 2. The noise correlator can by evaluated by ex- 
ploiting the 'flat interface' limit, valid to leading (zeroth) 
order in h. The result is 
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D 



(29) 



Equations (|28|) and ( p9| ) have precisely the forms antic- 
ipated earlier on general grounds. We note that, in the 
absence of thermal noise, our approach is very similar to 
that of Shinozaki Q . 

Finally, we have to impose the shear flow. To do this 
we write v = ^yex + u, where u is the deviation from 
the mean shear flow and should vanish far from the in- 



terface. Inserting this form for v in both ( |2 1|) and ( |22[ ), 
with the left side of ( p2|) set to zero appropriate to the 
viscous regime, we find that the shear term drops out of 
both the Navier-Stokes equation and the incompressibil- 
ity condition. We conclude that u plays exactly the same 
role in the sheared case as v plays in the unsheared case, 
and that the effect of the imposed shear is to add a term 
iyd x <j) to the left side of (p3|), just as in models A and 
B, and therefore a term (i/2)jk x [h 2 ]] l to the left side of 
(p8|), which then becomes 

-gf + ipk x [h ] k = —j— h k + Tfr{t) . 



(30) 



given by d c = (9 + /i)/2 for fi > 1, while for fi < 1 we 
will see that the situation is less clear. 
The scale transformation takes the form 

x = bx', x_L = 6 f x'_L, h = b x ti, t = b z t' . (33) 

To make further progress it is necessary to know whether 
£ < 1 or £ > 1. Since the shear term tends to enhance 
the interfacial coarsening in the x-direction, we expect to 
find C < 1 whenever the shear is relevant, though £ = 1 is 
possible for d > d c , where the shear rate 7 is formally an 
irrelevant variable. We will further argue that £ > 1 is 
unphysical, and will accordingly restrict consideration to 
C < 1 in the following. We will find, however, that the na- 
ture of the theory for d > d c differs according to whether 
fi > 1 or fj, < 1. We will therefore consider these two 
regimes separately. The former regime includes models 
A (/1 = I) and B (/1 = 2), while the latter includes model 
H (n — 0). A brief discussion, in the present context, of 
the case fj, = 1 can be found in pi). This special case 
had also been discussed earlier in the (physically very 



different) context of a sandpile model |19|. 



A. The case fj, > 1 



III. RENORMALIZATION GROUP ANALYSIS 

The starting point of the RG analysis is equation 
(|l|). Since, however, the system is anisotropic we ex- 
pect difference scaling properties in the directions paral- 
lel and perpendicular to the shear. Under coarse grain- 
ing, anisotropics will develop in the linear terms in the 
equation. Additionally, from the structure of the non- 
linear (shear) term it is clear that terms analytic in k 2 
will be generated in the response function self-energy and 
the renormalized noise. Anticipating this, we generalize 
equation ([!]) to (in Fourier space): 

d t h k + ^k x (h 2 ) k = -(A|k| 1+ ^ + u x k 2 x )h k + 7fc(t) . 

(31) 

The noise correlator takes the form 

(Vk(t)r)-k> (t')) = (Dlk^- 1 + D x k 2 x ) 5 k , k , 5(t - t') . (32) 

We apply a momentum-shell RG in which, for conve- 
nience, we impose an ultraviolet momentum cut-off, A, in 
the x-direction only. The RG transformation consists of 
three steps: (i) eliminating modes with A/6 < \k x \ < A 
(hard modes); (ii) rescaling the length scales, x and Xj_, 
the field variable, h, and the time, t; (iii) looking for fixed 
points of the equation of motion at which the theory is 
invariant under (i) and (ii). As usual, the elimination of 
modes will be executed perturbatively near the critical 
dimension, d c , of the theory. We will show that d c is 



A value of £ less than unity implies anisotropic scaling. 
Furthermore, in such cases the transverse part, k_L, of k 
dominates over k x in the terms involving powers of |k|, 
both in the equation of motion and the noise correlator, 
which then take the following forms: 

d t h k + \pk x {h 2 ) k = -(A|kJ 1+ ^ + v x k 2 x )h k + vu(t) . 

(34) 

(77k(i)?7-k'(0> = (Dlk^- 1 + D x k 2 x ) Sum' S(t - f) . 

(35) 

Note that for /i = 1 the term Ak 2 , coming from A|k| 2 can 
be absorbed into the v x k\ term, while the term Z?|kp I_1 
becomes a constant. So the case fi = 1 is covered by the 
general structure of equations (Em|) and (|35|). 

Applying the transformation ( p3| ) to equation (|^) then 
yields rescaled values for the parameters in the equation 
and the noise correlator: 

7 ' = 6* +2 - 1 7 (36) 

A' = b z - {1+ ^ A (37) 

v' x = b z - 2 v x + --- (38) 

B , = h z-2 X -i- {p -iK-(d-2K D (39) 

D' x = ^-2x-3-(d-2) C D x + ---, (40) 

where the ellipses indicate that the parameters v and 
D x acquire perturbative corrections due to the coarse- 
graining step of the RG procedure. By contrast, the pa- 
rameters 7, A and D acquire no perturbative corrections 
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- equations (p6|) , (37) and ( |39j ) are exact. The absence 
of perturbative corrections to 7 follows from the invari- 
ance of the general equation of motion, ([!]), under the 
transformation h — ► h + ho, x — > x + "{hot, which is the 
analog for our system of the usual Galilean invariance of 
Burgers equations (see, for example, pj). The absence 
of corrections to (|37]) and (|3^) follows from the fact that 
the vertex 7 carries a factor k x . As a result, all pertur- 
bative contributions to the response function self-energy 
and the noise correlator carry factors of k x . 

Let us first examine the linear theory (7 = 0) to iden- 
tify the critical dimension d c . In the linear theory, there 
are no perturbative corrections, and equations (p7|)~(|4C|) 
all hold exactly. From (|37|)-(|39|) we obtain 



zo = 2 , Co = 



1 + m 



7 - n - 2d 

xo = ^iTTir ' (41) 



where the subscripts indicate that these are the results of 
the free theory. Inserting these exponents into equation 
© gives D' x = b~ 4 ^ 1+f ^D x , indicating that D x flows 
to zero at this fixed point. 

Equation (^) determines the relevance, at the trivial 
fixed point, of the shear rate 7. From ( f4l| ) we obtain 
Xo + z — 1 = (9 +/x — 2d)/[2(l + fj,)] . Hence 7 is relevant 
for d < d r , where 



dc = (9 + A*)/2 , n>l 



(42) 



For d < d c , we expect a new fixed point to appears at 
which 7, A and D are all non-zero. Equations j3^), d37j ) 
and ( |39| ) give the corresponding exponents exactly: 



3(1+/*) 
6 + 2/i - d 



C 



6 + 2/i - d 



X 



3 — (j, — d 
6 + 2/i - d ' 

(43) 



We recall that in order for our calculation to be consistent 
we must have £ < 1, such that |k| ~ |kj_|. From relations 
(|l|) and © we see that this condition requires fi > 1, 
consistent with the case we are currently analyzing. 

Exponents ( fflf ) are correct only if the fixed point val- 
ues of the parameters 7, A and D are all non-zero, other- 
wise their scaling dimensions cannot be set equal to zero. 
To check this fact we perform a one-loop RG calculation 
to compute the perturbative corrections to v x . In gen- 
eral, integration over the hard modes gives the following 



equation for the renormalized propagator G < (k, u>) (see 
Fig.l): 

G< (k, u) = G(k, u) + G(k, u) E(k, u) G< (k, oj) , (44) 




FIG. 1. Dyson equation for the propagator in terms of the 
bare propagator (single lines) and the self-energy (hatched 
circle) . 



where the bare propagator is given by 

G(k,w)- 1 = -iu + v x k 2 x + A|k x | 1+M 



(45) 



and the self-energy S(k, oj) must be calculated perturba- 
tively in 7. From the relation, 



G < (k,oj)- 1 = G(k,w)- 1 - E(k,w) , 



(46) 



we clearly see that the perturbative corrections to v x 
come from terms of order k x in E(k, oj). Setting b = e , 



(47) 



with / infinitesimal, equations (|3q) and (Uq) yield, 



dv x 
~dJ 



(z-2)- lim 



f 



k^o v x k% I 



E(k,0) 



The standard one-loop diagram for the self-energy is 
shown in Figure 2 (see, for example, IhJ). 




q+k/2 




_(q+k/2) 



-q+k/2 



FIG. 2. One-loop contribution to the self-energy. The in- 
ternal lines also carry frequency labels (not shown). 

Full circles represent 7- vertices, open circles represent 
contractions of the noise, (fjk^-k) ~ D, and arrows are 
bare propagators. The leading term of the self-energy in 
the limit (k, u>) — > is given by 



E(k,0) = - 7 2 J D / 

Jn, q > 



1. 1 k i 



■q,n 



G(--q,0 



2(6 + 2/i -d) TT „ u2 



/i + 5 — d 



Uis x k 2 J (48) 



U 



Sd-2 



8(/i+ 1)(2tt) 



d-2 



6 + 2/i - d 



f 



q 3 — d — fj, — 

7 DX M+i i/j. 



6 + 
f + 1 



(49) 
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In the expression above T(u) is the gamma function 
and Sd—2 is the surface area of the unit sphere in d — 2 
dimensions. The notation (f2, q > ) means that we inte- 
grate with the measure c?Sl dq. x dq± / (2ir) d ~ 1 , within the 
outer shell Ae~ l < \q x \ < A. Due to the anisotropic na- 
ture of the non-linearity, there is no need to introduce 
a cut-off for q^. Furthermore, we have taken A = 1 
without loss of generality. 

Putting together equations fl47j), and and us- 
ing the scaling dimensions of the parameters 7, A and D, 
we finally obtain the RG flow equation for the effective 
coupling constant U, 



4 = 3 + 2/i, fi < 1 . 



(55) 



dU 9 + (i - 2d 



dl 



U 



2(6 + 2/i -d) 2 
(jjl + 1)(m + 5 -d) 



U 2 



(50) 



Consistent with our previous determination of the critical 
dimension, we see that the linear term in the flow equa- 
tion changes sign for d = d c = (9 + /i)/2. Moreover, the 
quadratic term is negative, implying that for any d < d c 
there is a non-zero stable fixed point U* = 0(e), with 
e = d c — d. The RG pcrturbative expansion is thus 
well behaved and the fixed point values of 7, A and D 
for d < d c are finite. The exponents (|4^) are therefore 
correct. On the other hand, for d > d c the only sta- 
ble fixed point is U* = 0, corresponding to an irrelevant 
non-linearity and thus giving the 'free' exponents of (41). 



Note that d' c differs from the critical dimension d c found 
for the case /t > 1 in (ff2|), namely d' c < d c . On the other 
hand they coincide in the limit [i — > 1. 

For d < d' c , from (^) one again obtains C < 1 an d 
therefore it is tempting to conclude that these are the 
correct exponents even for the /x < 1 case, provided that 
d < d' c . As a further consistency check, one may note 
that these exponents reproduce the ones of the free the- 
ory given by fl54| ) for d — > d' c . Unfortunately, the situation 
is not as simple as this. If we perform a one-loop per- 
turbative expansion below d' c , we formally get the same 
flow equation (|50|), since ( < 1 in this regime. How- 
ever, as we have seen, the fixed point of this equation is 
of order e = d c — d, which is not small for d ~ d' c . In 
other words, because of the gap between d' c and d c , the 
one-loop expansion in the form stated above is not under 
control in the regime /i < 1. We were not able to find 
a perturbatively consistent solution in this phase. As a 
consequence, we can only conjecture that the exponents 
we have found for /1 < 1 are correct, since they lack a 
substantial perturbative support. 

Finally, let us note that although one can formally find 
a solution with ( > 1, all the terms involving drop 
out at this fixed point, and the equation becomes essen- 
tially one-dimensional, which is unphysical. We therefore 
reject this possibility. 



B. The case /1 < 1 



For (j, < 1, equation ( [4l| ) gives C > 1 for the free the- 
ory, violating the assumption £ < 1 under which ( [fl| ) was 
derived. This suggests we look for a solution with £ > 1. 
In this case, k x will dominate over (or be the same order 
as) kj^ in |k|. The recursion relations for A, v x and D 
become 



A' = b z - (1+ ^ A 

u' x = b z - 2 v x + ■ ■ ■ 

D' = 5 z -2x-i-(M-i)-(d-2)c D 



(51) 
(52) 
(53) 



instead of (|37|) HJ39|) ■ At the fixed point of the free theory 
(7 = 0), equation ( |5l"| ) gives z = 1 + fi, so that ( |52| ) be- 
comes v' x = b^~ l v x , i.e. v x is driven to zero, since /1 < 1. 
The theory with v x = = 7 is completely isotropic, so 
Q = 1. Inserting z = 1 + ji and ( = 1 in (53) gives 
X = (3 — rf)/2. Summarising, the exponents of the free 
theory for /t < 1 are 



1 + M 1 Co 



Xo = (3 - d)/2 



(54) 



which coincides with ( p4| ) in the limit fj, — > 1. 

The relevance of 7 is again determined by equation 
(|36|). From (54), the combination \ + z — 1 is given, for 
the free theory, by \o + z o — 1 = (3 + 2/i — d) /2. Hence, 
for /1 < 1, 7 is relevant below the new critical dimension 



C. The case d = 2 

Some of the results derived above only hold for d > 2. 
This is because the idea that k_L dominates k x in |k| is 
clearly inapplicable in d = 2, since there is only k x . Sim- 
ilarly, the exponent £ can no longer be defined, so there 
are just two independent exponents, z and \. The equa- 
tion of motion and noise correlator are given by (|3~I| ) and 
(|35| ) respectively, but with |k| replaced by \k x \. We recall 
that model H (/1 = 0) is ill-defined for d = 2. 

The RG recursion relations for d = 2 become, 



7' = b x+z ~ 1 1 
A' = A 
v' x = b z - 2 v x + ■ ■ 
D' = b z - 2x -^ D 
D' x = b z - 2x ~ 3 D x 



(56) 
(57) 
(58) 
(59) 
(60) 



Equations (p6[), ( J5T[ ) and ( p9[ ) are exact, and therefore it 
seems that we have three equations for just two unknown 
exponents, x an d z. This apparent paradox is solved if 
one of the parameters is zero at the fixed point, since in 
this case the corresponding equation is trivially satisfied 
without setting the scaling dimension to zero, The shear 
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rate 7 is certainly relevant, since d = 2 is below the crit- 
ical dimension. If we assume D — > 0, using equations 
( |56| ) and ([57]) we get x + z = 1 and z = (J, + 1, giving 
X = — )U. This would imply a positive scaling dimension 
for D, which is inconsistent with D — > 0. Thus, we must 
assume A — ► 0, and find from equations ( ]56| ) and (59), 
X + z = 1 and z — 2% = /i at the fixed point, giving 



« = (2 + /i)/3, x=(l-M)/3 



(61) 



in = 2. Inserting these results into ( p7| ) gives A' — 
b~<~ 1+2 ^/ 3 \, so A flows to zero in d — 2, as assumed, for 
all n > -1/2. 



IV. STABILITY OF THE DOMAINS 

The calculations of the previous sections are important 
to assess the stability of the highly stretched domains in 
a coarsening system under shear. We recall that we are 
considering a shear velocity profile with flow in the x di- 
rection and gradient in the y direction. We denote by xj^ 
all the directions orthogonal to both x and y for d > 3. 
The effect of the shear is to stretch the coarsening do- 
mains, such that there are two different length scales, 
in, along the x direction, and L± in all the orthogonal 
directions. The transverse size of the domains, Lj_, is 
in general much smaller than longitudinal one, Ly |^],[|. 
What we have to check is whether the size A of the height 
fluctuation is larger than L± , inducing a breaking of the 
domains, or whether A < Lj_, meaning that the domains 
are stable under thermal fluctuations. 

In the long-time limit, the main orientation of the do- 
mains will be almost completely parallel to the shear flow, 
and therefore height fluctuations in the surface of the do- 
mains grow in a direction orthogonal to x. In d = 2, this 
implies that the only relevant fluctuations are in the y, 
that is h, direction. On the other hand, for d = 3, there 
are also fluctuations growing in the xi direction, which 
are not described by equation (Q). These two cases will 
therefore be treated separately. 



hn (i x/2 ,if) 



(63) 



In two dimensions we need only consider models A (/i 
1) and B(/i = 2). 



1. Model A 



In this case the critical dimension is d c = 5, so for 
d = 2 the shear is relevant. From the former sections we 
have x = an d z = 1. Equation ([33]) therefore implies 
that, whatever value L\\ takes, the height fluctuation A 
will be of order unity. In || it has been shown that 
for model A the transverse domain size is Lj_ ~ 0(1). 
This is an analytical result obtained in the context of the 
Ohta-Jasnow-Kawasaki approximation. This gives, 



A ~ Lj_ (d = 2 , Model A) 



(64) 



We conclude that model A in two dimensions is a 
marginal case, and we cannot exclude the possibility that 
thermal fluctuations may in this case break the domains, 
giving rise to a stationary state. 



2. Model B 



For model B we have d c — 11/5 > 2, and the exponents 
are x = —1/3, z — 4/3. Also in this case, therefore, we do 
not need to know the coarsening exponent for Ly, since 
from relation (|6|) it is clear that a negative value of x 
implies a saturation of A to a constant value: 



A -0(1) (d = 2 



Model B) 



(65) 



A. The case d — 2 

In two dimensions the height fluctuations of the sur- 
face are given by the fluctuations of the field h. Thus, as 
a consequence of the scaling relation h(x, t) = b x h'(x' , t') 
(see Eq. (|H|)), the height fluctuation A grows as 

A ~ h ~ i*/*F(t/Lf ) , (62) 

where the scaling function F goes to a constant for small 
argument and F(s) — s~ x ^ z for s — > 00. This means that 
if i 1 / 2 <C L\\ the surface grows like A ~ t x ^ z : whereas if 
t x l z we have A ~ L x . We can incorporate both 

limits in the form 



This result opens up two different scenarios, according to 
the the growth law for L±. If L± — i 1 / 3 , as argued in Q] 
by means of numerical experiments and RG arguments, 
then A <C L± and the domains must be stable against 
thermal fluctuations. If, however, L± — O(l), as sug- 
gested by some recent numerical simulations pi, then, 
as in model A, we cannot exclude the possibility that 
a breaking of the domains by thermal fluctuations oc- 
curs. Our result shows that a Lj_ — i 1 / 3 growth law and 
a thermally induced stretching and breaking mechanism 
are not compatible. Conversely, if a thermally-induced 
breaking of the domains is clearly observed in numeri- 
cal experiments, this strongly suggests that the relation 
L± - 0(1) holds. 
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B. The case d — 3 



1. Model A 



In three dimensions the situation is more complicated. 
First, as in d — 2, there are height fluctuations in the y 
direction, A y ~ h, described by equation (|l|). Secondly, 
there are fluctuations in the x± direction, Aj_, which can 
also become larger than L±, and that are not described 
by equation ([TJ). Thus, before assessing the stability of 
the domains for d = 3 we must formulate an equation for 
the description of these latter fluctuations. Fortunately, 
this will turn out to be a linear equation, such that no 
perturbative RG analysis is necessary. 

In order to describe surface fluctuations which grow in 
the x± direction we have to introduce a new height field 
h± which satisfies the equation, 



n 



(66) 



to be compared with (Q) . The operator C is still given at 
low momenta by C ~ A|k| 1+A \ Equation Jflq ) is linear, 
and therefore we can work out the exponents exactly by 
means of simple scaling. By setting, 

x = bx', y = b c y', h± = b x h' ± , t = b z t' , (67) 

and imposing scale invariance of equation ( |66] ) , we obtain 
(with the usual hypothesis £ < 1), 



A' = b z -^ +1 h 

D' = p-2x-C^-i_o 



(68) 
(69) 
(70) 



and setting to zero the scaling dimensions of all three 
parameters gives 



fi + l 
At + 2' 



C 



1 



/' 



A* + 1 
2(/j + 2) 



(71) 



Note that C is smaller than one, consistent with our as- 
sumption. We see that \ is negative for all the three 
interesting values of /i [p, = 0, 1, 2), meaning that height 
fluctuations along the xj_ direction are always finite, 
A x ~ 0(1). 

We have to assess now the physical importance of A y 
in the context of domain coarsening. From the usual 
scaling relations we get, 



A y ~ h~t*''F(t/Lfu t/L z / c ) 



(72) 



In general, evaluating the magnitude of A y from this rela- 
tion is quite subtle, as we need to compare the interfacial 
coarsening and equilibrium regimes in both the paral- 
lel and perpendicular directions. However, as we discuss 
below, in all cases of physical interest we have \ < 0, 
implying that the interfacial fluctuations saturate. 



In the case \x — 1 and d = 3, Eq. ([43]), with /i = 1, 
gives x — — 1 /5, and therefore A y ~ 0(1). For model A 
it was been found in [@] that L± ~ i 1 / 2 , giving 

(d = 3 , Model A) . (73) 



Ay < L j 



In model A, domains are therefore stable against thermal 
fluctuations. 



2. Model B 

In this case also the exponent x is negative: Eq. (|4 
with /i — 2 gives x = — 2/7, and z = 9/7, yielding 



0(1) (d = 3 



Model B) 



(74) 



Even though no analytical results or numerical simula- 
tions studies are available at the present time for model 
B in d = 3, we certainly expect L± to grow with time in 
this case, and therefore the domains to be stable. 



3. Model H 

For n < 1, as we have seen, we have a different critical 
dimension, given by Eq. (|55|), which is exactly three for 
(i = 0. This implies, using either ( |54| ) or (p3|), that x = 
and z = 1. Once again, this is the marginal case, with 



A y - 0(1) (d = 3, Model H) 



(75) 



V. SUMMARY 

Interfacial fluctuations have been investigated in sys- 
tems subjected to an external shear flow. Interfacial dy- 
namics appropriate to systems with non-conserved scalar 
order parameter ("model A"), conserved scalar order pa- 
rameter ("model B"), and conserved scalar order param- 
eter coupled to hydrodynamic flow ("model H") have 
been studied. In each case the interfacial dynamics is 
described by a similar equation, of the form ([!]) , where h 
is the local height of the interface and in which the eigen- 
value spectrum of the linear operator C has the form (m) . 
The models differ principally in the numerical value of 
the exponent /i, which is given by 1,2 and for models 
A, B and H respectively. 

The interface equations have the form of anisotropic 
noisy Burgers equations. In each case, exact renormal- 
ization group (RG) arguments determine the exponents 
z, C, and x that characterise the coarsening, anisotropy, 
and roughening of the interface respectively. In all cases, 
X < 0, implying that the thermally induced interfacial 
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width approaches a finite limit at infinite time. A con- 
sequence of this result is that the domain structure of 
a coarsening system under shear is stable against (suffi- 
ciently weak) thermal fluctuations. 

The general framework revealed by the exact RG rela- 
tions was supported by explicit one-loop calculations for 
H > 1. For n < 1, however, no one- loop equations con- 
sistent with the expected critical dimension, d! c = 3 + 2/z, 
could be derived. Whether this is just a technical diffi- 
culty, or signals some important physical difference be- 
tween the regimes \i > 1 and /i < 1 , merits further inves- 
tigation. 
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